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A note on the history (May 2006): In my book "Large Scale Dynamics of Interacting 
i/-) . Particles" [S] I refer to an unpublished note from early 1985 on the BBGKY hierarchy 

for hard spheres. My main point there was to provide a direct probabilistic proof 
for the time-integrated version of the hierarchy. Over recent years there has been 
repeated interest in this derivation, which encourages me to make my note public. 
I decided to leave it in its original form including likely inaccuracies. The work of 
R. Illner and M. Pulvirenti [IP] appeared in September 1985, see also the book by 



C. Cercignani, R. Illner, and M. Pulvirenti [CIP], who prove the same result using 



in 
o 

\Q [ special flow representation and methods from the theory of differential operators. 
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in Physics, Springer- Verlag, Heidelberg, 1991. 



[IP] R. Illner and M. Pulvirenti, A derivation of the BBGKY-hierarchy for hard 
sphere particle systems, Transport Theory and Stat. Phys. 16, 997-1012 (1987), 
preprint DM-388-IR, September 1985. 

[CIP] C. Cercignani, R. Illner, and M. Pulvirenti, The Mathematical Theory of 
Dilute Gases, Applied Mathematical Sciences 106, Springer- Verlag, New York, 1994. 

The aim of these notes is to establish the integrated form of the BBGKY hierarchy for 
hard spheres as used by O.E. Lanford pQ in his proof of the validity of the Boltzmann 
equation in the Boltzmann-Grad limit. The idea of a direct probabilistic proof is 
inspired by a paper of R. Lang and X.X. Nguyen 

We denote by Xj = {qj,Pj) 6 AxR 3 position and momentum of the j—th particle. 
The hard spheres have diameter a (and mass one). They are confined to the region 
A. A is bounded and has a "smooth" boundary OA. Conditions on <9A ensuring the 
existence of the hard sphere dynamics are given in the thesis of K. Alexander 
p. 13/14, and we assume the validity of these conditions here. We have exactly A" 
particles, j = 1, . . . , N. The n-particle phase space, n — 1, 2, . . . , N, is 

T„ = {(xi,...,x„) G (A x 1 3 )"| \qi-q\ > a/2 for all 

q G <9A, |ft - qj\ > a ,i, j = 1,. . . ,n,i ^ j} . (1) 



1 



In a collision of two hard spheres incoming and outgoing momenta transform 
into each other as 

Pi = Pi - &[u ■ fa - Pj)] , 

P'j = Pj + • (Pi ~ Pj)] » ( 2 ) 
i 7^ j, with Cj G S* 2 . At the wall particles are specularly reflected, 

P'j = Pj ~ 2 %i) [n{qj) ■ Pj] , (3) 

where n(qj) is the unit outward normal at the point of contact. For the construction 
of the hard sphere dynamics we refer to Alexander's thesis. 

We remove once and for all from T n the set of points which in the course of time 
run into either a grazing or a multiple collision. The phase space with these points 
removed is denoted by T* . T n \ T* has Lebesgue measure zero. Then, for all £ G R 
and for every point (x±, . . . , x n ) G T* , the flow 

t^T t {n \x u ...,x n )er* n (4) 

is well defined. In particular T* is invariant under Tj: n \ 

If incoming and outgoing momenta are identified (also at collisions with the 
wall), then is continuous in t, i.e. for all (x±, . . . ,x n ) G T* one has 

Jim T t (n) (x!, ...,x n ) = (xi, ...,x n ). (5) 

Here we will not identify incoming and outgoing momenta, i.e. we regard them as 
distinct phase points. For (xi, . . . ,x n ) G T* the map t i— > Tj n \xi, . . . ,x n ) is then 
piecewise continuous and we have to distinguish between the limit from the future 
(+) and from the past (— ) defined by 

T}± (x!,...,x n )= lim T^±l (xi,...,x n ). (6) 

e^0,e>0 

If the added signs ± are omitted, it is understood that the quantity in question is 
independent of the way the limit is taken. 

A function p n : T n — > R is continuous along trajectories of on T* , if for all 
[Xi, . . . , x n ) G T n 

lim p n (T t (n) (xi, ...,x n ))= p n (xi, ...,x n ), (7) 

where both the limit from the future and the past are understood. This implies then 
that for all (xi, . . . , Pi, . . . , (ft + acu, pj, . . . , x n ) G T* one has 

p n (x!, ...,q i ,p i ,...,q i + au,pj, ...,x n ) 

= p n (x u ...,q i ,p' i ,...,q i + au,p'j, ...,x n ) (8) 
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and similarly for collisions with the wall. 

For future convenience we define some sets: Let 

^N-n(Xl, • • • , x n ) = . . . , Xn) G rAT_ n | |<ft — qj \ > a, 

for i — 1, . . . , n and j — n + 1, . . . , iV} (9) 

for (xx, . . . , x n ) G T* and 

r N _ n (x 1 ,...,x n ) = Q (io) 

otherwise. Let 

flj(xx, . . .,x n ,p n+ x) = {uj e S 2 \(xi, . . .,x n ,qj + auj,p n+1 ) G C S 2 (11) 

for j = 1, . . . , n, (x\, . . . , x n ) G T* , and p n +i G M 3 and let 

Qjixx, . . .,x n ,p n+ x) = (12) 

otherwise. We define two subsets, fJj±, of Qj by 

Oj±(ari, • • . ,x n ,p n+ i) = {<2> G fij(xi, . . . , z n ,p n+ i)|u> • (p n+ i -pj) > (< 0)} . (13) 

After these preparations we can state our assumptions on the initial (t = 0) 
measure. 

Let P be the initial probability measure on Tat. P is assumed to satisfy: 

(i) P is symmetric in the particle labels. 

(ii) P has a density, 

P(dxi . . . dxiy) = /n(xi, . . . , xn)cLxi . . . dxjy- (14) 

(iii) /n is bounded by the canonical equilibrium distribution, i.e. there exist 
constants c, f3 > such that 

N 

f N (xi, ■ ■ ■ ,x N ) < c Y[hp(pj) (15) 

3=1 

on r^, where hp(p) = (^r) 3//2 e~^ p2//2 is the normalized Maxwellian. 

(iv) J'n = on Tn \ and J'n is continuous along trajectories of on Pj^, cf. 
© and ©. 

(v) The time evolved measure P t has a density /jv(£) given by = on r^T^ 
and 

f N (xi, ...,x N ,t) = f N (T^\xi, . . .,x N )) (16) 

for (xi, . . . , xn) G r^r. The canonical equilibrium measure is denoted by P oq . 

To avoid confusion we remark that identities are always understood pointwise. 
If they hold only a.s., we state this explicitly. Often we will work with densities of 
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measures. As in the case of /jv(t) and of p n if) below we will choose then a specific 
version. 

Because of a definite number of particles, N, the correlation functions are, up to 
a multiplicative factor, just the marginal measures. We fix a particular version of 
these measures by 

p n (x u . . . , x n , t) = N . . . (N - n + 1) / dx n+ i . . . dx n f N (x 1 , ...,x N ,t). 

Jr N _ n (xi,...,x n ) 

(17) 

Note that p n {i) = on T n \T* by our definition of T N _ n (xi, ■ ■ ■ , x n)- Let A C T„ be 
a Borel set. We remove a set of Lebesgue measure zero to guarantee that A C T* . 
Then because of the hard core exclusion and by symmetry 

/ dx\ . . . dx n p n (x\, • • • , x n , t) 
J A 

= N...(N-n+ l)P{( Xl (t), . . . , x n (t)) e A} , (18) 
where for j — 1, . . . , N we set 

x j {t±,x) = {Tg ) x) j . (19) 

The probability in ()18|) is independent of whether the limit is taken from the future 
of from the past. 

To avoid an overburdened language it is convenient to set 

t > 0, 

which we do from now on. This is no restriction, of course. 

Let r m > 0, m = 1, 2, . . . , be the time of the m-th collision between the set 
of particles with labels 1, . . . , n and the set of particles with labels n + 1, . . . , N. If 
there are simultaneous collisions between the two groups of particles, then they are 
ordered according to the label in the first group. 

Proposition 1 The following idendity holds for all Borel sets A C T* n , for all 
n = l,2,...,iV, 

/ dx\ . . . dx n p n (x\, • • • , x n , t) 
J A 

dxi . . . dx n p n (Tzj/ fa, ...,x n )) 

oo 

+ N ■ ■ ■ ( N - n + !) p { Mr m +), • • • , x n (r m +)) G T^_ t+ A, r m < t} 

m=l 

CO 

- ^ JV . . . (N - n + l)P{ (x!(r m -), . . . , x n {r m -)) G T^_ t _A, r m <t}. (20) 

m=l 
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Proof: We use inclusion - exclusion to obtain 

P{( Xl (t),...,x n (t)) e A} 
= P{(x 1 (t-),...,x n (t-)) e A,n <t} 
+P{(x 1 {t+),...,x n {t+)) e A,n >t} 

oo 

= p {( x i(t-)i-- -,x n {t-)) G A,r m+ i > t,r m < t) 

m=l 

+P{(x 1 ,...,x n ) eT ( y_A,r 1 >t} 

oo 

= P{(x 1 (r m +), . . .,x n {r m +)) G T^_ t+ A,r m+1 >t,r m < t} 

m=l 

+P{(x!, . . . , x n ) e t5 } A} - p{( Xl , . . . , x n ) e Tit A, n < t} 

oo 

= P{(Mr m +), ■ ■ .,x n (r m +)) G rW_ t+ A,r m < t) 



m=l 

oo 



-J2 P{(xi(r m +),...,x n (r m +)) e^_ i+ A,r m+1 < i} (*) 

m=l 

+p{(x 1 ,...,x n ) erl?A} 

oo 

- ^ P{(x!(r m -), . . . ,x n (r m -)) G 7^_ t _A, r m < t} 

m=l 

oo 

+ ]T P{(xx(r m+ i— ), . . . , x n (r m+1 — )) G P^ +1 _ t _A, r m+1 < t} . (**)(21) 



m=l 



To justify ()21|) we need an integrable bound. Clearly, the m-th term is bounded by 
cP eq {r m < t}. We will show in Lemma El below that this bound is summable. 
In (J2D) (*) and (**) cancel each other because as sets 

{(xi, . . .,x N ) G T* N \(x 1 (T m +), . . . ,x n (r m +)) G T^_ t+ A, r m+1 < t} 
= {(a? x , . . . ,x N ) G T* N \(x 1 (T m ^i-), ■ ■ .,x n (r m+1 -)) G T^ +i _ t _A,r m+1 < t} .(22) 

□ 

Lemma 2 Let r m > 0, m = 1, 2, . . . , be the time of the m-th collision for the system 
of N hard spheres (collisions with the wall are not counted). Then 

oo 

P eq{r m < t} 

t / dqidpx \ dp 2 / dQj a 2 uj ■ (p 1 - p 2 )Pc q ,2(<?i,£>i, <?i + au,p 2 ) ■ (23) 



m=l 
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Proof: We think of the hard sphere dynamics as a flow under a function (special 
flow), cf. [3] for this construction in our context, and we prove Lemma |2] first for 
this case. 

Let B be the base and T : B — ► B be an invertible map which preserves the 
finite measure fi. Let h : B — > M + be the ceiling function. We assume that h is 
integrable. The phase space is then T = {x <E B,y E M+ |0 < y < h(x)}. The 
flow T t is constructed piecewise in the following way: T t : (x, y) i— > (x, y + t) until 
the first time for which y + t = h(x). Then (x, /i(x)) i— ► (Tx,0). We refer to this 
transformation as a collision. The construction is then continued into the future 
and the past. The measure fi(dx) x dy = P eq is invariant under T t . Let r m > 0, 
m = 1, 2, . . . , be the time of the m-th collision. Then we claim that 

oo 

P ^{ T ™ ^ = *M#) • ( 24 ) 

m=l 

Since h > 0, Yl'jLo h(T~ix) = oo fi(dx) a.s. by the Poincare recurrence theorem. 
Therefore 

k-2 fe-l 

£ fc = {x G B| h(T- j x) <t,J2 h(T~ j x) > t}, (25) 

3=0 j=0 

k — 1, 2, ... , forms a partition of 5. Then 

oo oo „ 



Y P {- T m <t} = t/i(B 1 ) + Y v(dx)h(c 

m=l fc=2 

oo „ m—2 oo „ 

+ £{/ Kdx)(t-J2 h (T- j x))+ J2 / M^)MT- fe+2 x)} 

m=2 ^ Bm j=0 k=m+l J Bk 

oo 

= tJ2v( B k)=tK B )- (26) 
fc=i 

For hard spheres the base consists of configurations with outgoing momenta and 
is defined by 

B = {(xi, . . . ,xn) € T* N \ there exists a pair i ^ j, 

such that qj = qi + aou, (pj — pi) ■ lj > 0}. (27) 

The ceiling function is defined as the time until the next collision (not counting 
collisions with the wall) . The equilibrium measure induces on B the invariant surface 
measure 

N N N 

{ Y a 2 dqiduj (b ■ (pj - p { ) ] [ dq k ) JJ hp(p k )dp k . (28) 

i^=j=l k=l,kj^i,j k=l 

Its total weight is given by (|23|). □ 

We want to express (|2U|) in terms of correlation functions. For this purpose we 
first have to show some regularity of these functions. 
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Lemma 3 Under our assumptions on P, for every s£l, p n {s) =0 onT n \T* n and 
p n (s) is continuous along trajectories ofT^ on T*. 

Proof: Since, by assumptions (iv) and (v), /at(s) has the same continuity properties 
as /jv(0), we may set s = 0. 

To simplify notation we abbreviate x — (xi, . . . , x n ),y = (x n+ i, . . . ,xn) and we 
set x(t±, x) = Tj:±x. For x G T* let A(x, t) C A be spatial region traced out by the 
particles' motion x(s),0 < s < t, with initial condition x. We set A(x, 0) = A(x). 
Correspondingly we define A(y,t) C A for y G Then for a; G T* let 

T w _ n (x,t) = G IVJAOM) n A(y,t) = 0}. (29) 

We have T N _ n (x) = T N _ n (x, 0) up to a set of cfo/-measure zero. For x G T* we define 
the flow on rN_ n (:r)* of N — n particles in the spatial region A \ A(x). Here 
the * indicates again that we remove from T N _ n (x) a set of Lebesgue measure zero 
on which the flow remains undefined. 
With these definitions, for x G T*, 



PniT^x) = [ d yN...(N-n + l)f N (T t ^x, y) 

Jr N _ n {x(t)) 

[ dyN...(N-n+ l)f N (T t ^x, T^y) (30) 

JT N . n (x{t)) 



by Liouville's theorem for the map T t {x{t)) for fixed t. 

We choose now a r such that < t < r. For x G T* and y G Fjv_„(x, r) the 
"x" -particles and the "y" -particles do not interact during the time interval [0, r]. 
We then have two possibilities: 

(1) The time evolution exists into the future and the past, i.e. (x,y) G T* N . In this 
case (T^x^T^^y) = T^\x,y). We denote the set of such y's by T N ^ n (x, r). 

(2) The time evolution does not exist, i.e. (x,y) T* N . In this case, by assumption 
(v), 

f N (T t (n) x,T t {x(t)) y) = for < t < r. (31) 
Therefore, for every x G T*, 

\Pn(T^x) - p n (x)\ 

= N...(N-n + l) [ dyf N (T t m (x,y)) 

dyf N {T^xM xm y)- [ dyf N (x,y) 



<N...(N-n + l){ dy\f N (Tr(x,y))-f N (x,y)\ 
+c dyf eq , N (x,y) + c dyf e<ljN (x, 

v) • (32) 

Jr N - n (x(r))\r N - n (x,T) Jr N - n (x)\r N - n (x,T) 



The last two terms are bounded by const, r. For fixed r the first term vanishes in 
the limit t — ► 0. This follows from dominated convergence and our assumption (iv). 

□ 

Lemma 4 Under our assumptions on P, for every (xi, . . . , x n ) G T* the map t \— > 
p n {xi, ■ ■ ■ ,x n ,t) is continuous, i.e. 

lim p n (x!, . . . ,x n ,s + t) = p n (xi, .. .,x n ,s) . (33) 

Proof: We use the same notation as in the proof of LemmaE! Since, by assumptions 
(iv) and (v), Jn(s) has the same continuity properties as /W(0), we may set s = 0. 
For every x G T* we have 

p n {x,t) - p n {x) 



dyfN(x,y,t)- dyf N (x,y) 
{y\(x,y)er* N } J {y\(x,y)er* N } 

dy(f N (T^\x,y))-f N (x,y)). (34) 

{y\(x,y)&r* N } 

The claim follows then by dominated convergence from assumption (iv). □ 

Proposition 5 The following identity holds for every Borel set A C T* n , n = 
1,2,. ..,N, 

dx\ . . . dx n p n {x\, . . . , x n , t) 

dx\ . . . dx n p n (T_ t {x\, . . . , x n )) 



A 

n nt r 

+ ^ / ds dxi...dx n [C jin+ ip n+1 (s)](T^ +s (x 1 ,...,x n )). (35) 

j=1 Jo J A 

Here the collision operator is defined by 

(C j ^ n+ ip n+1 (s))(x 1 ,...,x n ) (36) 
= a 2 / dp n+1 / duu ■ (p n+1 - Pj)p n+ i(xi, ...,x n , qj + au,p n+1 , s) . 

J JQ(xi,...,x n ,p n+1 ) 

Proof: We consider the third term of ({ZD}) , cf. Proposition d For < s < t we 
want to compute the limit as e — > of 

^2 -P{(xi(T m -), . . . ,x n (r m -)) G T™_ t _A,r m G [s, s + e]} 

X 

~P S {particles with labels 1, . . . , n collide exactly k times with 



m=l 



fc =l £ 



particles with labels n + 1, . . . , N during the time interval [0, e], 

at the times r of collision (x\(t— ), . . . , x n (r— )) G T^ r _ f _A} . (37) 
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Here P s is the measure P evolved to time s. By assumption (iii) the sum for k > 2 
is bounded by 



oo 

C 



k-P eq {particles have exactly A; collisions during the time interval [0,e]}. (38) 



e 

k=2 



By the same argument as in Lemma El 
lim - P cq {particles have exactly one collision during the time interval [0, e]} 

= / dq x dp x \ dp 2 / dQj a 2 u ■ (p l - p 2 )Peq,2(ft,Pi, ft + au,p 2 ) (39) 

and according to Lemma |21 

oo 

kP cq {particles have exactly k collisions during the time interval [0, e]} 

k=l 

= e dqidpi \ dp 2 / du a 2 uj ■ (p x - ^2)^,2(91, Pi, ft + au,p 2 ) ■ (40) 

J J J0l-(X1,P2) 

Therefore in the limit e — > the expression in ()38|) vanishes. 

We are left with the term k = 1 of (p?7J). Let us label the particle at the collision 
with n + 1. Then we have to compute the limit e — > of 

(iV — n)-P s { particles with label 1, . . . , n collide exactly once with particle n + 1 
and do not collide with particles with label n + 2, . . . , N during the time 

interval [0, e], at the time r of collision (a?i(r— ), . . . ,x n (r— )) 6 Tf"^_ t _A} 
n x 

= ^^(^V — n)-P s { the only collision during time interval [0,e] is 
3=1 ^ 

between particle j and particle n + 1, at time r of collision 

Mr-), . . . , x n (r-)) e T^ T _ t _A} + 0(e) . (41) 

The error is bounded by 

-P eq {there is more than one collision during the time interval [0, e]| , (42) 

which vanishes in the limit e — ► 0. 
Let 

A^e) = {(x l ,...,x n+1 )eT* n+1 \T^ n+l \x 1 ,...,x n+1 ) for 0<t<e 

has as only collision the one between particles j and n + 1, 

at time r of collision {x\(t— ), . . . , x„(r— )) e T s ^ T _ t _A}. (43) 



In the definition flUJ) Xj(t±, x) = (T t i ± +1) x) j . If the set defined in (gT]) is called B^e) 
then 

n 1 - 

(gH) = VViV - n)- / dxi . . . dx N f N (xi, ...,x N ,s) 
j= i e Jb 3 {s) 

n 

= ]>>-n)- \ 

+ >(iV — n)- / dx% . . . dxNf]sr(xi, . . . ,xn, s) 

j= l el JB 3 (e) 

- dxi...dx n+ i dx n+2 ...dx N f N (x u . . .,x N ,s) . (44) 

JAofe) ir»_„_i(ii,...,i„ + i) 



dxi . . . dx n+1 / dx n+2 ■ ■ ■ dx N f N (xi, x N , s) 

Aj(e) JV N - n -2(xi,...,x n+ i) 



The second term is again bounded by (J42|) and vanishes therefore in the limit e — » 0. 
Multiplying with the factor N . . . (N — n + 1) of (|20|) we are left with 



n 1 f 

/] - I dxi . . . dx n+ ip n+1 (xi, . . . ,x n+1 ,s) . 
, =1 e 



(45) 



Let r, < r < e, be the time of collision. Then on Aj(e) 



q n +i = qj + au + r(pj - p n+1 ) . (46) 
We perform this substitution in the integral ()45j) . The change in volume element is 
dqjdpjdq n+1 dp n+1 = a 2 uj ■ (pj - p n+l )dq 3 dpjdTdQdp n+ i . (47) 
We flow on Aj(e) all n + 1 coordinates from time to time r. Then 

n 

2_\ ~ I dx x ... dx n+1 p n+1 (x 1 , . . . , x n+1 , s) 

j=l £ J Me) 

= - / dra 2 I dx\ . . . dx n l dp n+ i I duu ■ (pj — p n+ i) 

7= 1 £ JO JT { s ll_ t _A J JQ :j .(x 1 ,...,x n ,p n+1 ) 



3 

,(n+l) 



><Xh(^i, • • -,x n ,uJ,p n+ i)p n+ i(Tr T >(x u . . . } x n} qj + au,p n+1 ),s) 
- / dra 2 I dx\ . . . dx n f dp n+ i f duu ■ (pj — p n+ i) 

7= 1 E JO JT { s ll_ t _A J Jn^(x 1 ,...,x n ,p n+1 ) 



3 

,(n+l) 



xp n+1 (T\ > (x 1} ... } x n ,q j + au3,p n+1 ),s) + 0(s) . (48) 

In the second integral xs is the indicator function of the set {xi, . . . , x n , u, p n +i\ 
T^ +1 \xi, . . . , x n , qj + au,p n+ i) for — r < r' < e — r has only one collision}. [As 
collisions we always refer to collisions between two particles and not to collisions 
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with the wall. Therefore in ()46)1 we should actually use the free flow of particles 
j and n + 1 separately including collisions with the wall. After flowing the n + 1 
coordinates to time r we still obtain (|48)).] The error term is again bounded by (J42j) . 
Only P e q refers now to the equilibrium measure of n + 1 particles. 

To obtain the limit as e —>■ of (|4*Hj) we have to show that the integrand is 
continuous at r = 0. To see this we bound as 



dx 1 ... dx n / dp n+1 / du;u ■ (pj - p n+1 ) 

™' T _ t A J Jn j -(x 1 ,...,x n ,p n+1 ) 

xp n+ i(Tl" +1) (xi, ...,x n ,q j + au,p n+ i),s) 

/ dxx . . . dx n \ dp n+ i / dCjuj ■ (pj - p n+x ) 

JtHa J Jn i .(x 1 ,...,x n , Pn+1 ) 



'«_t«-» " JClj- (xi,...,X„,Pn+l) 

xp n +i{xi, ...,x n , qj + au,p n+1 , s) 



< / dxi . . . dx n I dp n+ \ / duo 

xw ■ (pj -p n+ i)p n+ i(T^ +1) (xi, . . .,x n ,qj + au,p n+1 ), s) 

+ / dxx... dx n j dp n+ i / duu ■ (pj - p n+ i) (49) 

Jt { s ™\& j Jn j _(x 1 ,...,x n , Pn+1 ) 

x|p n+1 (Tl" +1) (xi, . . . ,x n , qj + auj,p n+ i), s) - p n +i(xi, ■ ■ • ,£n,<?j + aw,p n+1 , s) \ . 

(s) 

In the first term we bound p n + x by const. / e q,n+i- By dominated convergence this term 
vanishes then in the limit r — > 0. In the second term we integrate only over points 
such that (x\, . . . , x n , + auj,p n+ i) G Therefore by LemmaElthe integrand is 

continuous in r and vanishes as r — > 0. 

Altogether we have shown that the measure 

oo 

^ JV . . . (N - n + l)P{{ Xl {r m -), x n (r m -)) G T^_ t _A, r m G ds} (50) 

m=l 

is absolutely continuous with respect to the Lebesgue measure and has a density 
given by 

a 2 / dx x ... dx n / dp n+1 / du>u> • (p, - p n+1 ) 

j=i J T s-\ A J Jn j .(x 1 ,...,x n , Pn+1 ) 

xp n+1 ( Xl , ...,x n , qj + (£j,p n+1 , s) . (51) 

We note that by Lemma 0] this density is continuous in s. 

The same argument applied to the second term of (|2(J|) shows that 

oo 

J2 N . . . (N - n + l)P{(a; 1 (r m +), . . . , x n (r m +)) G T^_ t+ A, r m G ds} (52) 

m=l 
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has a density given by 

^2 a 2 dxx . . . dx n / dp n+1 / dQjQj ■ (pj - p n+l ) 

j = l J {(x 1 ,...,q j ,p' j ,...,x n )£T < s ™\A} J Jnj_(xi,...,x n ,p n+ i) 

xp n +i{xi, ■ ■ ■ ,x n , qj + au,p n+1 , s) 



= 2_. a 2 / dxi . . . dqjdp'j . . . dx n / dp' n+1 (53) 

j=1 J{( Xl ,...,q j ,p' j ,...,x n )eT£\A} ' ' J 

x / duu-(p' n+1 -p' j )p n+1 (x 1 ,...,q j ,p' j ,...,x n ,q j + au,p' n+1 ,s), 

where we used again Lemma H3 which ensures that on the domain of integration 
p^li is continuous through a collision. We relabel in (|53~j) (p'j,p' n+ i) as (pj,p n +i) an d 
subtract (jSH) from (JHSD - Since A C r* , 

n pt p 

y~] ds dx 1 ...dx n [C jn+1 p n+ i(s)](x l ,...,x n ) 

n pt p 

= ds I dx i--- dx n [C jtn+1 p n+ i(s)}(T^\(xi, ...,x n )). (54) 

Jo J A 



□ 



To obtain the integrated from of the BBGKY hierarchy we have to iterate (|35j). 
For this purpose we go back to (jUJ). Since we integrate there over a Borel set of r n+1 , 
we could have chosen any other version of p n+ \(s), i.e. any other function p„ + i(s) 
such that p n+ i(s) = p n+1 (s) dx\ . . . dx n+ i a.s.. We used however certain properties of 
Pn+i(s) in the proof below (j43|) . Therefore, if we want to replace p n+ i(s) by p n+ i(s), 
the latter has to satisfy: 

(1) p n+ i{s) = p n+ i(s) a.s.. 

(2) For fixed s, p n+ i(s) is continuous along trajectories of T^ n+1 ' on 

(3) For every (x lf . . . , x n+1 ) G r* +1 , s i-> p n+ i{x u . . . , x n+1 , s) is continuous. 

(4) There exist constants c', f3 such that 

p n+ i(s)<c'f^ n+1 . (55) 
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Corollary 6 Let p n+ i(s) : T n+ i — > R satisfy the Properties (QJ) to O) gwen above. 
Then 

I dx\ . . . dx n p n (x\^ . . . , x n , t) 
Ja 

= / da? 1 ...dx n /? n (Tl" ) (a;i,...,a; n )) 
7a 

+ ds dxi... dx n [C j>n+1 p n+1 (s)](T^ +8 (xi, 

j=Q JO J A 

Lemma 7 Let p n (t) be defined by 

Pnip'L . • • • j X n ^ t) p n (T_^ {x\i ■ ■ • i ^n)) 
™ /-t 



,£„)). (56) 



+ / ^[C iin+ ip n+ i(s)](T 

.7=1 ^ 



-t+s+(^l' • • • ) -^n)) (57) 



for every point (xi, . . . ,x n ) G T*, where p n+ i( s ) satisfies the above Properties (EJ) to 
Q). T/ien p n (t) satisfies also the above four properties. 

We note that p n {t) = on T n \ T* by the definition of p n and of C^+i. 

Proof: Property (JIJ follows from Corollary |3 Property (j3J) follows from Lemma El 
(continuity of p n along trajectories of on T* ) and from Property (jU) of p n+ i(s). 
For Property (0J) we note that by Assumption (iii) on the initial measure 

Pn(t)<c'fW n . (58 ) 

Together with the assumed bound on p„ +1 (s) this gives a bound of the desired form 
with some new constants c",(3'. 

For Property (0) the first term of (J5T|) is continuous along trajectories by Lemma 
121 Therefore we only have to consider the second term. We have 

\-t+( X 1i • • • i x n)) 



V/ ds[C^ n+l p n+l (s)]{T^l 
' i=i y ° 

rfs [Cj,n+lPn+l(s)](T5^ + (xi,...,2; n )) 

<^2\{ ds+ ds}\[C j!n+1 p n+1 (s-r)](T^ t ) +s+ (xi,...,x n ))\\ 
j=i Jo Jt 

iPn+i(s - 7") - C jin+ ip n+ i(s)](T_ t+s+ (xi, . . . ,x„))| . (59) 
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By assumption the integrand of the first term is bounded by c"/cq,n- This implies 
that the first term vanishes in the limit r — > 0. In the second term we use dom- 
inated convergence. By the assumed continuity of s i— * p n+ %(xi, . . . , x n+ i, s) for 
) G r* +1 the integrand vanishes pointwise in the limit r — > 0. □ 

By Lemma [7| we may set in (jHBj) p„ +1 (s) = p„+i( s ) and obtain 
dx l ... dxnPnixx, ...,x n ,t)= / dx 1 . . . dx n p n (T {n t ' + (x 1 , . . . , x n )) (60) 

A J A L 

n nt 

+ V / dt 1 (C jun+1 {p n+1 o T^MT^^+fa, x n )) 
h=i Jo 

n n+1 » t rt x 

+ dtl dt 2 {C jun+1 (C h ,n+2Pn+2 

We iterate N — n times and obtain 

Proposition 8 The following identity holds for every Borel set A C T* n , n = 

1 v. 

dxi...dx n p n (xi,...,x n ,t)= (61) 

A 

N—n n n+m-1 „t rt m -i r 

^2 ^2 ' ' ' S dt i - - dtm dxx . . . dx n 

m=0n=l jm =l Jo Jo 

X (Cji,n+1 • • • {Cj m . n +m{Pn+m ° ?t m + )) ° ' ' (^i -i+ ( Xl ' ' ' ' ' ' 

In a way Proposition |H1 is our final result. There are however two reasons for 
reorganizing somewhat the integral (|61|). First of all we would like to get rid of the 
continuity assumptions, i.e. we would like to extend the validity of (}6T| to a more 
general class of initial measures. Secondly the Boltzmann-Grad limit is not quite 
apparent in the given form of (|61jl. 

We introduce the notion of a collision history. I choose this name to distinguish 
it from a sequence of real collisions of the iV-particle system. One should remember 
that the correlation functions are averaged quantities. The correspondence between 
collision histories and sequences of real collisions is only very indirect. 

A collision history is specified by the following list: 

(a) n E N, 

(b) m e NU {0}, 

[If, as assumed so far, the number of particles equals N, then 1 < n < N and 
< m < N -7i.] 

(c) (x u . . . ,x n ) e r*, 
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(d) O er 

with the constraint < £ m < . . . < t± < t, m > 1, 

(e) Uu-.-Jm) eN m 

with the constraint 1 < ji < n, . . . , 1 < j m < n + m — 1, m > 1, 

(f) (p!,...,p m )eR 3m , m>l, 

(g) (i lr .,^)e(S 2 r,m>l, 

with a complicated constraint depending on n,m,Xi, . . . , x n , ti . . . , t m , ji, ■ ■ ■ , j m , 
pi,...,p m which is defined below. For future convenience we abbreviate a colli- 
sion history as (xi, . . -,x n ,5), where 5 stands for (m,ti, . . . , t m ,ji, . . . , j m ,Pi, ■ ■ -,Pm, 

Ld 1 , . . . , Ld m ). 

Given the collision history (x 1 , . . . ,x n ,5) we construct an evolution of particles 
in the following way. We choose n particles at (x±, . . . , x n ) G T*. We consider 
this as the phase point at time t, (x±, . . . ,x n ) = (xi(t), . . . ,x n (t)) and we evolve 
backwards in time up to t = 0. We evolve the phase point (xi(t), . . . ,x n (t)) to 
T_t +tl+ (x\(t), . . . , x n (t)) = (xi(ti) : . . . , x n (ti)). We add a particle with label n + 1 at 
q^iti) + aui with momentum pi. We require that a)i G Qj^xifa), . . . ,x n (ti),p n+ i), 
i.e. (xi(ti), . . . , x n (ti), q^iti) + acu,pi) G We call this new phase point of 

n + 1 particles (xi(ti), . . . ,x n+ i(ti)) and evolve it to T ( ^H 2+ (xi(ti) , . . . , x n+ i(ti)) = 
(xi(t 2 ), ■ ■ ■ , x n+ i(t 2 )). We add a particle with label n + 2 at q^^) + acj 2 with 
momentum p 2 . We require that Cj 2 G Qj 2 (x 1 (t 2 ), . . . , x n+ i(t 2 ),p n+2 ). We call this 
new phase point of n + 2 particles (xi(t 2 ), . . . , x n+2 {t 2 )). The final step is to evolve 
(^i (^m) ; • • • ; x n+m (t rn )) G r* +m to 

T -tt+ > ( x ^( t rn), x n+m (t m )) = (xi(0), . . . , x n+m (0)). If m = 0, then we only evolve 
( Xl (t),...,x n (t)) = ( Xl ,...,x n ) toT^ m \ Xl {t),...,x n {t)) = ( Xl (0),...,x n (0)). To 
make the dependence on xi, . . . , x n , 5 explicit we write Xk(s,xi, . . . , x n , 5) for < 
s < t, in particular Xk(xi, . . . , x n , 5) = Xfc(0), k = 1, . . . , n + m. 
Let 

A(xi, ...,s n ; [0,t]) 

be the space of all collision histories for given n, starting configuration (x 1 , . . . , x n ) G 
T* and time span [0, t}. A(xi, ...,s n ; [0, £]) is a subset of 

J |J ... (J (I X I 3 x S 2 )" 1 
m>0ji=l im=l 

defined by the above construction. We define a measure dS on A(xi, . . . , x n ; [0, t]) : 
c?(5 is the counting measure with respect to the discrete indices m,ji,... ,j m and the 
Lebesgue measure otherwise. 
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Given a collision history (x\, . . . , x n , 6) we define the weight function by 

TO 

W{x u ...,x n ,5) = YlW^k • (Pk~ Pj k (t k ,x 1 , . . .,x n ,5))} . (62) 



k=l 



Lemma 9 The following identity holds for every Borel set A C T* n , n G N, 

/ dx\ . . . dx n p n (xi, . . . , x n , t) = dx\ . . . dx n I d5 (63) 

JA JA JA(xi,...,x n ;[0,t]) 

xW(xi, . . . ,x n ,5)p n+m ( S ) (xi(xi, ...,x n ,5),.. . ,x n+m ( S )(xi, . . .,x n ,5)) . 

Proof: We write out (|5T|) using the definition of Cj >n+ %. From Assumption (iii) it 
follows that 

Pn+m — (' feq,n+m ■ (64) 

(The simpler part of) Lanford's estimate on the uniform, in a 2 N, convergence of the 
BBGKY hierarchy shows that 



dx\...dx n \ d5\W(xi, . . . , x n , S)\ 

A J A(x,...,x n ;[0,t]) 

X feq!n+m(8)( Xl ( Xl i ■ ■ ■ ■> X ni $\ ■ ■ ■ i X n +m(S) {Xl , • • • , X n , 5)) 

< c'" / dxi... dx n f^'l( x i, ■■■,%n) (65) 

with f3' < (3 for all t. The details of this estimate can be found in F. King's thesis 
jS]. Therefore the integrations in (JBTj) may be interchanged freely. □ 

In the form (jo^j) we can extend our identity to a more general class of initial 
measures. In particular we will remove the restriction of a definite number of par- 
ticles. 

Let T be the grand canonical phase space, 

r = |J r n . (66) 

n>0 

The grand canonical equilibrium measure with inverse temperature (3 > and fu- 
gacity z > is defined by 

1 1 " 1 

f!*' t n( x U ■ ■ ■ > x n )—dxi ...dx n = — \\{zh p (p j )}—dxi ...dx n (67) 
n\ Z - LJ - n\ 

i=i 

on T n ,n = 0, 1, ... , where Z is the normalization constant. Let C be the class of 
functions / : T — > M such that 

(i) f n is measurable, 
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(ii) f n is symmetric in the particle labels, 

(iii) there exist positive constants M, z, (3 such that 

\fn(xi, ...,x n )\< Mf^\ Xl , ...,x n ) (68) 

for all (xi,...,x n ) G r n , n = 0, 1, . . .. 

Note that actually T n = for sufficiently large n because of the hard core 
exclusion. 

Given /eCwe define the "correlation function vector" p : T — > R by 
°° 1 f 

Pni%li ■ • ■ i *£n) / , T / d^n+l ■ ■ ■ d^n+mfn+mi^li ■ ■ ■ > •^n+m) (69) 

^— ' m! / r f Tl r n 
for all (xi, . . . , x n ) G T n , n = 0, 1, 

Lemma 10 Let T : C — > C be the map defined by \69) . Then T is one-to-one and 
onto. 

Proof: We have by assumption (iii) 

IPra^l; • • • 5 %n)\ — M ^ ^ j" / dXn+i . . . dx n _|_ m / e q n _|_ m (Xi, . . . , X n + m ) 
m=0 ' J 

<Me^f^( Xl ,...,x n ). (70) 
The inverse map is given by 

f ( \ (-^ m f 
f n (x u ...,x n ) = 2^—^r / 1 



dx n +i . . . cfe n _|_ m p n _|_ m (xi, . . . , x n _|_ m ) . (71) 



m=0 



□ 

Let P be a signed measure on T with density /. Then the time evolved measure 
P t has the density fit) with one version given by 

f n (xi, ...,X n ,t) = f n (T^ + (xi, X n )) (72) 

for all (x\, . . . , x n ) G T* and f n [t) = on T n \ T* , n — 1, 2, — If / G C then also 
/(t) G C and the correlation functions at time t are still defined by ()69|) . 

Theorem 11 Lei P be a signed measure on T with density f G C. Then for every 
Borel set A C Y n , n — 1,2, . . ., the following identity holds 

dx\ . . . dx n p n {x\, . . . , x n , t) 



dx\...dx n I d8W(xi, . . . ,x n ,5) 

J A(xi,...,x n ;[0,t]) 

X Pn+m{S){x\{xi, . . . , X n , 5), . . . (xx,...,x n ,5)). (73) 

The right (left) hand side of does not depend on the chosen version of p [pit)). 
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Corollary 12 Under the same assumptions, independently of the chosen versions 
of p{t) and p, 



Pn ipbl j • • • j t)= dSW(xi, ...,x n ,S) 

JA(xi,...,Xn,[0,t]) 

XPn+m(S)(Xi(x 1 , ...,X n ,8),.. . , X n+m ( S ) (x x , . . -,X n ,8)) (74) 

dx\ . . . dx n a.s.. 

Proof: Since the sum over m is finite, we consider a term with fixed m for reasons 
of notational simplicity. We abbreviate x — (xi, . . . ,x n ). Then ()73j) reads 



dxp n ix,t) = dx d5W(x,5)p n+m (y(x,S)) (75) 

J A* JA(x;[0,t],m) 

with the obvious definition of y(x, 5) G r* +m . Here we replaced A by A* = A n T* 
which leaves the integral unchanged. Let C* C C be the class of densities / which 
satisfy the continuity assumptions (iv) and (v), cf. (jIBJ), omitting the requirements 
of normalization, positivity, and definite number of particles. If / e C*, then (J75|) 
holds by Lemma 01 For an arbitrary / e C there exists a sequence f £ G C such that 



lim/ £ = / a.s. (76) 



on T. Consequently 
and, since 
we also have 
Therefore 



lim p £ = p (77) 



e->0 



lira f it) = fit) a.s., (76 



limp e (t) = pit) a.s.. (79) 



£^0 



lim / dxp e n ix,t) = / dxp n (x,t) . (80) 
J a J A 

We have to investigate now the convergence of the right hand side of ([75)1 . Let 
p G C* and p = a.s.. Then also p(t) = a.s. and we conclude that 



0=/ dx d5Wix,5)p n+m iy(x,S)). (81) 

JA* J A(x;[0,t],m) 

Since all sets of measure zero can be approximated in this way, we conclude that 
y^ 1 (considered as a mapping for sets) maps sets in T* n+m of dx x . . . ofe n+m -Hieasure 
zero to sets in {(x, 6)\x G A*, 5 G A(x; [0, t], m)} of dxcM-measure zero. 
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By (|ZZj) there exists a set T C T* +m such that T n+m \ T has measure zero and 
such that limp^ +m = p n+m pointwise on T. Let Xf t> e the indicator function of the 

set T. Then 

Pn+mXf 

(y(x,6)) (82) 

and by the argument given above 

Pn+mi 1 ~ Xf)(y( X ' 5)) = dxdS a ' S - • (83) 

Together with (jfiHj) our claim follows from dominated convergence. □ 
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